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Faddeev equations in configuration space and integral form for three-atom scattering processes 
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is displayed. This formulation appears to be a valuable alternative to current approaches based on 
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I. INTRODUCTION 

The Faddeev scheme for three bodies is a very successful framework in nuclear physics and specifically in the 3- 
nucleon system 0,0- Most of the calculations are performed in momentum space, which is naturally adapted to 
nuclear forces derived from meson exchange diagrams. Also the coordinate space version of the Faddev equations in 
differential form has been successfully used, especially for the three-nucleon bound states Q , but an integral equation 
form of the Faddeev equation in configuration space has not yet been developed in detail. To the best of our knowledge 
this has not been done before, neither in the fields of nuclear nor atomic physics. 

While in nuclear physics there is only one 2-body bound state, the deuteron, diatomic molecules have numerous 
vibrational and rotational levels. This poses quite a numerical challenge, especially in 3-atom scattering processes. A 
great deal of work has been done in this field using hyperspherical harmonic expansion methods, |ld|.|ll|. Here and 
below we quote only a few recent papers. It is the purpose of this paper to describe an alternative approach using 
the configuration space Faddeev equations in integral form. We have chosen that scheme since over the years a highly 
efficient numerical technique in coordinate space has been developed by one of us |J] for solving 2-body scattering and 
bound state problems. This method is based on spectral expansion techniques in terms of Chebyshev polynomialsi 
denoted as S-IEM in Q. As we shall indicate these methods are ideally suited to the integral Faddeev scheme for 3 
atoms. It is hoped that this technique will result in an accuracy for the final 3-body results that is higher than what 
is currently achieved in nuclear physics calculations in momentum space ( about 3 significant figures). Needless to 
say there exist a great demand for theoretical support for various kinds of 3-atom processes: 3-atom recombination 
processes 12], collisional cooling of co-trapped atomic and molecular ions by ultracold atoms[l4] | .[T ^| .[T8l|.[T^|. quantum 
dynamics of ultra cold atom-diatom collisions |l5j . Three-atom potential surfaces have been developed, for instance 
for H3\ij% or for Li3 [l7|. 

In section II we provide a brief derivation of the form of Faddeev equations we want to use. The coordinate space 
representation expressed in terms of vectors is subsequently set up in section III, while its partial wave representation 
is worked out in section IV. Because of the complexity of the resulting expressions for general orbital angular momenta 
a separate section V is devoted to an overall s-wave reduction, since this case provides a more transparent insight into 
the structure of the equations. Technical details for the partial wave decomposition are deferred to the Appendices. 
We delegate a numerical feasibility study to a forthcoming paper, but we nevertheless add some remarks on that issue 
in section VI. Finally in section VII we provide a summary and conclusions. 



II. THE THEORETICAL FRAMEWORK 

Let us regard a system of three distinguishable atoms, which interact by 2-and 3- atom forces. Though apparently 
most often the full potential surface for 3 atoms is used ab initio, we prefer to separate out the 2-atom forces from the 
full potential surface and treat the remainder as a genuine 3-atom force. Thus we denote the total potential energy 
between 3 atoms by 

V = Vi + V 2 + V 3 + V 4 (1) 

where the 2- atom forces are conveniently denoted by Vi = Vjk (with j, k ^ i) and V4 is the genuine 3- atom force. 
Instead of introducing four Faddeev components corresponding to the four contributions to the potential we split 
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the 3- atom force into 3 parts, which are then combined with the three 2-atom forces, as follows. We introduce the 
notation 

Vi = V} 1] + Vf } + V} 3) (2) 

(i) 

where V£ , i = 1,2,3 has the property of beeing composed of two pieces, one symmetrical and the other antisymmet- 
rical in the interchange of atoms j,k ^ i. In detail 

V 4 = K(l,2,3) = ^[F(123) + F(132)] + i[F(123)-F(132)] (3) 
+ ^[F(123) + V(321)] + J[F(123) - V(321)] 
+ ^[V(123) + 1/(213)] + J[F(123) - V(213)] 

where the argument 123 stands for the space coordinate vectors of the 3 atoms. Thus the two terms in the first row 
are symmetrical and antisymmetrical, respectively, under exchange of atoms 2 and 3 and similarly for the two terms 
in next two rows. We group the potential into three terms Vi and in the 3-atom Schrodinger equation (SE), 
which then can be written as 

3 

(H -E)y = -Y / (V 1 + Vf ) )*. (4) 

i=i 

In the present paper we describe a scattering process where one atom is incident on a molecular bound state of 
the other two atoms. We can number the three atoms such that the target molecule is formed out of atoms 2 and 3. 
Thus the initial channel state is 

|01 >= \u >23 Wo >i (5) 

composed of the 2-atom bound state \u > with quantum numbers to be specified later and a momentum eigenstate 
Wo >i °f relative motion of atom 1 and the molecule (2,3). That channel state obeys 

(H Q + V 1 )\cf> 1 >=E\ ( t> 1 > (6) 

where E is the total center of mass energy, and Hq is the kinetic energy operator for the three atoms. 
The integral form of the SE is 

* = G (^(y 4 + n w )*) (7) 

where where Go is the free 3-atom propagator, Go = (E + ie — Hq)^ 1 . There is no extra driving term, since there 
is no solution to the left hand side of Eq.QJ alone which is regular at the origin and purely outgoing. The above 
equation suggests the decomposition 

3 

* = y>< (8) 



with 



^; = Go(^+F 4 W )* (9) 

Inserting the decomposition |JSJ into the right hand side of Eq. © and using the well known identities 

(1 - GoV^GoVi = G h (10) 
(l-Go^i)^ 1 ' =1 + G ii (11) 
following standard steps || one ends up with the 3 coupled Faddeev equations(FE) 

= 0! + Go ti(ih + tfo) + (1 + G h)G V^^ 
V>2 = Go * 2 (V>3 + V>i) + (1 + G t 2 )G V} 2) V 

i>3 = Go * 3 (V»i + tfa) + (1 + G t 3 )G V} 3) V (12) 
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The new ingredient, the 2- atom t-operators ij embedded into the 3-body space, with i = 1,2, 3, obey the Lippmann 
Schwinger equation(LSE) 

U = Vi + ViGoti (13) 

This generalization is natural, since the Faddeev equations contain a two-atom potential Vi separately in each 
arrangement i. A short description of the related two-body operator r(r, r') and the amplitude T (r) is as follows: 

The Lippmann Schwinger equation for the solution of the two-body SE is tp = <f) + 9oVip, where go is the Green's 
function (E 2 +ie — Hq2)~ 1 , Hoi is the two-body kinetic energy operator, V the potential, and <fi is the incident plane 
wave. The product V(r) ip(r) is denoted by T(r), and is expressed in terms of an integral involving the r operator as 

/>oo 

V{r)t/}(r)= T(r,r') 4>(r') dr'r' 2 . (14) 



By manipulating operator identities (I1L)|) and Ijlll) one can show that (1 — Vgo)~ 1 V = r, which is equivalent to Eq. 
(|13|l . with the difference that go contains only the kinetic energy for the relative motion of two particles, while Go 
also includes the kinetic energy for the other particle, and further, E 2 is the two-body energy, while E is the total 
three-body energy. The connection between and the two-body t— matrix r is established in section III, see Eqs 156(1 
and J57J). 

The kernel pieces in Ecis. <|12|) determine the asymptotic behavior of the wave function components ipi in the outgoing 
channels. For the 3- atom break- up they can be read off from the expressions standing to the right of the free 
propagator Go 0] ■ 

Thus if we define the three T-amplitudes as 

ipi = <fii+G Ti 
-02 = G T 2 

03 = Go T 3 (15) 

then we obtain for the Tj 's the expressions 

T\ = iiG (T 2 + T 3 ) + (1 + tiG )G F 4 (1) * (16) 

and correspondingly for T 2 and T3. The T- amplitudes are more useful than the wave function amplitudes because 
from the T- amplitudes one can derive the asymptotic normalization of the wave functions and furthermore , most 
importantly, the T's decrease to zero at large distances. Inserting Eas. <|15H into the right hand sides of Ea. Hltj|) and 
the ones for T 2 and T3 one finds a set of 3 coupled equations 

T x = hG (T 2 + T 3 ) + (1 + tjGo) VP [0! + Go (T x +T 2 + T 3 )] (17) 

T 2 = Mi + (l+taGo) n (2) [^i +hG (T 3 +T 1 )] 

+ (l + t 2 Go)F 4 (2) Go(T 1 +r 2 +r 3 ) (18) 

T 3 = t 3 <h + (1 + t 3 G ) V} 3) [fa + i 3 G (Ti + T 2 )] 

+ (l + t 3 G )F 4 (3) G (T 1 +T 2 +T 3 ) (19) 

This is the set of equations which are to be solved. A partial wave decomposition will be given in section IV. 
Because of (jHJ the total break- up amplitude will be 

U = T 1 +T 2 + T 3 (20) 

from which the physical break- up amplitude can be obtained by calculating the matrix element 

<0o\U o > (21) 

where \<fia > is a product of momentum eigenstates for the free motion of the three atoms in the final state. 

One is also interested in the transition amplitudes for elastic and rearrangement scattering. They can be read off 
from Eas. (|12l> by rewriting the kernel parts such that Vi distorts the Green operators 

G t = (22) 

E + iE-H -Vi y ' 
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Since G$ti = GiVi , with (1 + GqU)Gq — Gi, one finds for the elastic transition amplitude 

U n = W 2 + i> 3 ) + v} 1) «f 

and for the two rearrangement amplitudes 

U21 = V 2 (^ 3 + ^i) + n (2) ^ 

u 31 = v 3 {ip 1 + i> 2 ) + v} 3) v 

Upon inserting (|15f) into the above expressions for the U's, we obtain 

C/11 - VlG {T 2 + T 3 ) + y 4 (1) [0x + G (Ti + T 2 + T 3 )] 

C/21 = ^[0i + G (T 3 + Ti)] + U 4 (2) [0! + G (T! + T 2 + T 3 )} 



(23) 

(24) 
(25) 

(26) 
(27) 



Uai = V 3 [<f>i + G (Ti + T a )] + U 4 (3) [<£i + Go(2i + T 2 + 21 



The physical amplitudes are 



where the 0j's obey the SE 



< 4>i\Un > 



Vi\<Pi > = (E-H )\</>i > 



(28) 
(29) 
(30) 



Like in Eq. (J5J the channel states 2 and cj> 3 contain bound states. Thus when calculating the matrix elements l|29l) 
one can replace Vi<f>i by Gg *0i, and hence replace ViGo in Eas. l(26T) . ljT7)l . (j2*5)l by 1. These substitutions then lead to 
the final expressions 



Un = (T 2 + T 3 ) + V} 1} [0i + Go (7\ + T 2 + T 3 )] 



(31) 



C/21 = Gq Vi + (T 3 + Ti) + U 4 (2) [0! + Go (Ti + T 2 + T 3 )} 



(32) 



U 31 = Go Vi + m + la)] + ^i 3) [<£i + Go (T x + T 2 + T ; 



(33) 



It can easily be verified that the forms in the equations above for the Un are identical to the more standard expressions, 
which are 



Ua = (V j + V k + V i )% j,k^i 



For instance one has 



< 0i|Z7n > = < 0i 

= < 0i 

= < 01 

= < 01 

+ < 01 

= < 01 



V^ 2 + ^ 3 ) > + <0i|V 4 (1) |*> 

G c 7 1 (^2+^ 3 )> + <0l|U 4 (1) |*> 

Gq 1 !* > - < 0i|GoVi > + < 0iln (1) l* > 

^1 + ^2 + ^3 + F 4 |* > - < 0i|G o " 1 G o (T4 + U 4 (1) )|* > 

U 4 (1) |*> 

V 2 + U 3 + V 4 (1) |*> 



(34) 
(35) 



In the second to last step we used the definition of the Faddeev amplitudes, Eq. 

In the case of identical atoms the 3 Faddeev amplitudes © are identical in form, only the particles are permuted 
8] . One easily finds 



"02 + 1p3 = 'Pi'l 



(36) 
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where 

V =Pi 2j P 23 + P13P21 (37) 
is a sum of a cyclical and anticyclical permutations of 3 objects. Thus the total state can be written as 

$=(1 + V)fa (38) 

and hence only one FE is needed that reads 

fa = 0i + GohVfa + (1 + Go<i)G ^ 4 (1) (1 + V)fa (39) 

? If one defines again 

fa = fa + G T (40) 

then the amplitude T obeys 

T = tiPfa + (l + t 1 G )vj; 1) (l + P)fa 

+ t 1 VG T+(l + t 1 G )v} 1 \l+V)G T (41) 

The driving term in 1)41(1 is contained in the first line and the integral operator is contained in the second line. The 
complete break up amplitude is 

U = {1+V)T (42) 

Because of the identity of the atoms there is only one amplitude for the transition into two-body fragmentation, 
which, according to Ea. ((39[) is 

U = ViVfa+VP(l + V)fa 

= VG Q 1 fa+VT + vj 1 \l + V)fa + vj 1 \l+V)G T (43) 

The last expression is valid for calculating the physical matrix element < fa \ U > . 

This concludes the derivation of the formal framework. The case where only two atoms arc identical is similar to 
the above and is left to the reader. 

III. COORDINATE SPACE REPRESENTATION 

In a 3-atom system there are three 2-body fragmentation ( or arrangement ) channels going with three types of 
Jacobi vectors ( ijk = 123 and cyclical permutations) 

-» -* -{i) -* ™j -* (AA\ 

X — Xj Xki y — %i ' Xj ■ Xk V J 

rrij + rrik rrij + mu 

We introduce coordinate space states 

\x y >!= IfW,^ > (45) 
\x y >2= \x (2 \y< 2) > (46) 
|f y> 3 = |f (3) ,y< 3) > (47) 

Each set of states is complete: 

J dxdy \xy >< xy\ = 1. (48) 

The various terms for the T-amplitudes in Eqs. (|17fl - ( I19|l will now be written in the coordinate space represen- 
tation, using the definitions in Eqs. I|45() - ( I47|) . It is natural to represent the amplitude T\ as 1 < xy\T\ >, T 2 as 
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2 < xy\T2 > and T3 as 3 < xy\T 3 >. Then after inserting completeness relations the first term on the right hand side 
of Ea. lfr7|) becomes 

1 < x y\hG \T 2 + T 3 > = J dxdy < x y\h\xy >i J dx'dy'i < x y\G \x"y" >i (49) 

J dx"'dy"'[i< x" y"W" r > 2 2 < 3f" y"'\T 2 > 
+ x < x" f\x>" f> > 3 3 < x'" f'\T 3 >] 

Various matrix elements occur which we now consider one by one. The 2-body t-matrix obeys the Lippmann Schwinger 
equation which in configuration space reads 

< x y\t\x" y 1 > = < x y\V\x" if > (50) 
+ J dx"dy" < x y\V\x" y" > J dx"'df < x" f\G Q \x"' y"' X x"' y"'\t\x' y 1 > 

The two-atom force V is diagonal in the vector y which means 

< xylVlx'y 1 >= 6(y- y 1 ) < x\V\x" >= -^-^ J dqe^-^ < x\V\x" > (51) 

The three-body free Green's function can be related to a two-body free Green's function, by making use of the property 
that each is diagonal in momentum space. We use Jacobi momenta related to x^ l \ y^ 1 ^ and obtain 

< x y\G \x' g > = I dp dq< x y\p q> -5- < p q\x' y 1 > 

/ P7 _i_ 7 v 2 2 — 

1 f dp dq e ip<s-s') e iq-{y-v') 1 



-1-3 J dqe^y-y^g{x,x'-e q ) (52) 



where g(x,x';e q ) is the well known free single particle Greens function 



x ; e g ) = p — — (53) 

2n \x — x'\ 



We introduced the reduced mass /ii of the two atoms in the arrangement (i) 



Mi = ^ (54) 

and Mi the reduced mass between the particle i and the pair ( jk) 

Mi = (55) 

and e 9 = E — jjj- is the energy related to the 2- atom subsystem. 

Similarily, we can relate the two-body t-matrix embedded in the three-body space to the two-body r-matrix defined 
in the two-body space as follows. Inserting (|5T|) and 1(521) into the LSE lj50)l it follows that the t-matrix has the form 

< xyltlx'y' >= j^ys J dqe^y-y^ < x\r(e q )\x' > (56) 

where t(x,x'; e q ) =< x\T(e q )\x' > obeys the two-body LSE 

T(x,x';e q ) = V{x)5(x~ x*)+ [ dx"V(x)g{x, x"; e g )r(x", x'; e q ), (57) 
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and where we have assumed that V(x) is a local potential. 

The recoupling matrix elements in (|49ll requires some consideration. The three different sets of Jacobi vectors are 
linearily related to each other. Thus x^ and ifi 2 ^ in arrangement 2 can be expressed in terms of x^ and y^ 1 ^ in 
arrangement 1 as 



Similarity 



where 



x^ = Ax w + By<^ (58) 
y< 2 ) = Cx (1) + (59) 

f 0) = A > s (i) + B itfi) (60) 
y<3) = C'x^ + D'tf^ (61) 

A = ^— B = -l (62) 

m 2 + m 3 

C = m 3 m i + m 2 + "t 3 ^ _ mi 

7712 + 7773 TOl + 771 3 777-1 + 777-3 

A' = ^ — B' = l (64) 

"72 + 777 3 

C l _ 777 2 7771 + 777 2 + ^3 jj, _ ™l ^ 

m2 + TU3 mi + W2 rni + 7772 

As a consequence one can express the state \x y >i as 

\x y>i= \Ax + By, Cx + Dy > 2 (66) 
\x y > x = \A'x + B'y,C'x + D'y> 3 (67) 

The above equations mean that a spatial configuration represented in arrangement 1 by vectors x and y is represented 
in arrangement 2 by vectors Ax + By and Cx + Dy and similarily for arrangement 3. This leads to 

1 < x y\x' y 1 > 2 = 6{x' -Ax- By)8{y^ -Cx- Dy) (68) 
1 < x y\x' ff > 3 = 8{x' - A'x - B'yjS^ - Cx - D'y) (69) 

Using now EqsgHl,©, ©, EHJ, EHJ we obtain 

1 < x y\h Go \T 2 + T 3 >= ^3 J dqj dy" J dx" (70) 

[ 2 < Af" + By", Cx" + Dy"\T 2 > +3 < A'x" + B'f, Cx" + D'y"\T 3 >} 

Thus the unknown amplitudes 2 < x y\T2 > and 3 < x y\T 3 > occur with shifted arguments under the integral. 
Next we regard the second term on the right hand side of Eq. (|17l) : 



1 <xy\(l+ h G )U 4 (1 Vi >= J dx'dg 1 <xy\vj 1) \x'y' > xx <3> &\<h. > (71) 
+ J dx'dy 1 , < x g[hG Q )\tf #>! J dx"dy" 1 < x' y'\v} 1) \x" y" > u < x" y"\<j> x > x 

The 3- atom force V4 has to be a scalar under rotations. Therefore for spinless atoms and assuming locality it 
will have the form 

1 <x y\V^\x' y 1 > 1 ^V^{x 1 y,x ■ y)5{x - x')5{y - y 1 ) (72) 
Therefore, together with lf72 |) .l(55 |) .lJ52 )l we obtain 

<x y\(l +t 1 G Q )v} 1) |0i > =VW(i l y ) i.i)u(S)^ I e < M 

+ J2WI d( lJ dS ' I dtfe*V-W J dx"r(x, x"; e q ) 
xg(x",x";e g )VU(x',y',x' . y'^^j^e^ (73) 
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Finally the third term on the right hand side of Eq. (|17|) can be transformed by means of similar steps with the 
result 

<xy\(l + tiG ) V} 1) G (T 1 +T 2 +T 3 )> (74) 
= V^\x,y,x-y)-^ [ dqj dx' [ dy'e i ™-^g(xx';e q ) 



w 

(< x'y'ln > + <Ax' + By 1 , Cx' + D$\T 2 > 
+ < A'x' + B'tf,C'x' + D'g\T 3 >) 

+ 7^3 J dqj djf'e*™-*"* J dx" J dx'r(x,x , ;e q )g(x , ,x";e q ) 

V W(x",y",x"-y")-^- s J dtfj dx'" J dfe^'^'h^"^) 

(< f "f'\T x > + < Ax"' + By"', Cx'" + Dy""\T 2 > 
+ < A'x:'" + B'y"', Cx'" + D'y"'\T 3 >) 

We leave it to the reader to work out the coordinate space representations for the remaining two equations 1)18(1 . ( 

us). 

In case of identical particles there is only one FE given by (|41ll to be solved for the amplitude T(x, y) =< xy\T >. 
This eqation is of the form 

T = T° + K (75) 
where T° is the driving term and K incorporates the T- amplitudes. In the coordinate space representation 

T\x,y) = (1)3(^)9/2 | J d fj d ^Hy-y) (76) 

[r(x, \f + e q )u no {- 2 -y" ~ -j?) + r(x, ~f - j^; e>„ (j^' + \f)]e^" 
+ VW (x, y, x ■ y) [u no (x)e^ 



and 



+ u no (-^x-y)e^<i s -^ + ^(-itf+jfle'ft-C-f*-**)] 

+ T^Tg f dc lf d y" J dx"e lq ' {y ~ y ) J dxT(x,x;e q )g(x,x";e q ) 

VU(x",y",x" ■ r)j^[(u no {x")e^y' 
+ u no {-\x" - y")e^<i s "-^ + u no {-\x" + f) e ^<-i s "-^')] 

K(xy) = J dx" J dy" J dqe^~^ J da?r(x , a? ; e q )g(a? , a> '; e q ) (77) 

[T{-tS" y", U" tf) + T(- V + y", -±3» - ±f)] 

+ V ^(x,y,x-y)j^ J dqj djfe*&-*) J dx'g(x,x';e q )[< x'$\T > 

+ <Ax' + By 1 , Cx' + Dy\T > + < A'x' + B'y', Cx' + D'$\T >] 

+ J dqj df J dx'e^y-y^ J dx"T(x,x";e q )g(x'',x';e q )V (1 \x',y',x' -y') 

^ dqe tq ~y ^ g{^' ^" ' 



,x ,t q i^ 



(2tt) e 

[< x"f\T > + < Ax" + By", Cx*' + Dy"\T > 
<A'x" + B , y",C'x" + D'y"\T>] 
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Once the amplitude T(x, y) has been determined the physical 3— atom break up matrix element is obtained as 

<pq\U Q >= <pq](l + V)T> (78) 
' dx dy e -*^* e -*f ? T(x, y) + permuted parts (79) 



(2tt) 3 

The permuted parts are most conveniently evaluated by applying the permutation V to the left: 

i <pq\V=2<P q\+3 <pcfr (80) 

Then one expresses the Jacobi momenta of the type 1 in terms of the Jacobi momenta of the types 2 and 3. One has 
in case of identical particlesjg 

^) = -i^) + 3^ 2) (81) 

^ = -p< 2 ) - ~g< 2) (82) 

p<D = _^3) _ ^3) (83) 

g<l) = p{3) _ 1 -<3) (84) 



and therefore 
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2 <#?1 =i< -^P+ ( 85 ) 
and similarily for 3 < pg| . Thus altogether 

< p glC/o >= T(j?fl + T(-ip + - A q, -p- ±q) + T{-\p- -^p- X -q) (86) 

It remains to display the physical matrix element < 4>'i\U > for elastic and inelastic scattering , where \fa >= \u' > 
|<f > represents the final state. According to (|43|l one obtains 

< fa\U > = < fa\PG^\fa > + < 4>' X \V\T > (87) 
+ < fa\vP{l+V)\fa > + < fa\\V^(l+V)G T> 

For the first term we can use the SE for fa and obtain 

< falVG^lfa > = < fa'lVVlfa > (88) 

= J dx'dy'dx"dy" < fa^x'tf >i i < x*y'\V\x"y" >u < x"y"\V\fa > 

This is an example where one might evaluate the V -matrix element differently from the one shown above. The 
expressions <|68|) . <|69ll can be rewritten as 

1 <xy\vWy'>i = (^) 3 s(x - y - g^ff - -if - AD ~ BC v) m 



, 1 x8 .,. 1 ^ D' s .^, A' _ A'D'-B'C 
,-^) 3 S(x - -y 1 - -y)6(x< - -g , 



For identical atoms one has A = -\,B = —1,C = |,D = -\,A' = -\,B' = 1,C = —j,D' = -\. n Then we 
obtain 

< fa\VG^\fa > = < fa\VV\fa > (90) 



' K^) 3 / dy I dge-^-Ve^-y' (91) 



(2tt) 3V 3^ 
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This together with the remaining terms on the right hand side of 187(1 yields 

< &\U >= 7^(^) 3 / d vJ dge^-ye^'-y (92) 

K'i^y 1 + §0)^(1 + \y\)u na {-\$ - ~y) + - |s)v-(||^ + \y\WA\y 1 + \y)\ 

+ ^e-^^vO^i- y)K (x)e^ 

+ u no (-^x-y)e^<i s -^ + u no {~x + 0) e **>-(-i*-i*)] 

(< fy" |T > + < + By", Cf + Dif\T > + < A'x' + B'tf, C'x' + D'$\T >) 
For distinguishable particles corresponding results occur and their evaluation is left to the reader. 



IV. PARTIAL WAVE DECOMPOSITION 

In view of the various vibration-rotational levels it is natural to decompose the integral equations for the various 
T-amplitudes into partial waves. We exemplifiy that step for the case of three indentical atoms, and choose therefore 
Eas. l|75|l . (|76|) . H77|l . The scalar quantities r as well as the single particle Green's function <?, both contained in Eqs. 
(|75[) . (|76|) . l|77[) have the partial wave decomposition 

T(x,x';e q ) = ^2 Yi, m (x)Ti(x,x';e q ) Y£ m {x' ) (93) 
g(x, x 1 ; e q ) = ^ Y hm (x)gi(x, x';e q ) Yj*„(£') (94) 

I, in 

and 

Stf-i?) = fc^l Y? <m {x') (95) 

xx f — ' 

Here v is the unit vector pointing into the direction of v. When used as an argument of a spherical harmonics it stands 
for the angles 8 and <fi. 

As a consequence the LSE for r, Ea. (|57|l . for a given angular momentum I reads 

n(x,x';e q ) = V(x) 5{x ~ f ] 
xx' 

+ V(x) dx"x" 2 9l (x,x";e q ) Ti(x",x';e q ) (96) 
Jo 

Further ingredients are the initial bound state 

(*) Y htmo (x) (97) 

and the standard expansion of the plane wave 

e* = 4vr i X jx (qy)Yl, (q)Y x ^(y). (98) 

Next one introduces bi-polar spherical harmonics of total angular momentum L 

y^ M (x, y) = J2 °( lXL > m»M)Y ltm (x) Y x ^{y). (99) 

m,fi 
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where the C(- • • ) is a Clebsch- Gordan coefficients as defined in p|. However, in what follows we are going to omit 
the third magnetic quantum number, since m + fi = M . 

Using the notation above we will now begin to evaluate the amplitude T given in Eqs. 1)41(1 and 1(75(1 . ((76(1 . ((77() . The 
two driving terms in the first line of Ea ((41|) will be denoted by T 0,1 and T ' 2 , respectively. 

As an intermediate result one obtains for T 0,1 

< xylT ' 1 > = < x y\tV\<t> > (100) 



"3' (2ir) 9/ 

\ An ) l\ LM 

x |°° dqq 2 j x (qy) J ' d y> J df n ( q y')y^ M ^y" + ^y',y / ) 

4 2 2 4 > 

*n(x,\-y' + -y\;e q ) u no (\-y' + -y\) E i Xo jx {Qov")Yx 0<flo [qo) 

E C(?oAoL',mo,Mo)^ A ? +A10 (^y" + [R'° + (-)'], (101) 

n Next we expand n(x, \ W' + £g)u no (\ly" + into Legendre polynomials, using P k (x-y) = {-) k ^=yl '°(x, y) 

V k ' 

and obtain 

n(x, \\f + y\;e q )u no {\\f + = ]T 2nV2k + T(-) k y° k '° k (y',y")G k , (102) 



where 

Gk = J dt P k (t) 7j(ac,, 3 .., .,.•,■,■-.,.■•„,,.,._,;,■ >- 
and where t is the cosine of the angle between the vectors y " and y ' After placing qo into the z direction we obtain 



G k = ! + dtP k {t)T l {x^' + y\-,e q )u no {\lf + \y*\), (103) 



<xy\tV\<fi> = (1)3^^(4^^^ y?; x M {x,y) 

(27T' ' 



IX LM 



X/ dqq 2 j x {qy) dy'iy' 2 d y" 2 Jx (qy') 



X fdy' [ dy"y; L x M {\y" + \y',y') E (^)VWTT(-) k y° k ° k (y\y'')G k 



x 



3" 3' 
E 2 Ao JAo(WV(2Ao + l)/(47r) 



x E C(/ AoL',mo,0) y^™ ^?' + [(-)'" + (")']■ (104) 



L' 

The remaining angular integrations can be performed analytically with the result 



dy' dy''y; L x M {-y'' + -y' 1 y') W (£', y")^ + ^j/', 0") 



= $l,l< 6 M , mo rf£ loXo (y',y") (105) 

where the explicit expression for h is given in Appendix A. 

The r- matrix occuring in Eg ((103(1 contains as a driving term the function V (x) 8(x — x')/x x' ( see Eg. (19611 ). It 
is more convenient to explicitly separate this term out by defining the function r; 

n {x, x'; e q ) = V (x) felil + ^(a;, x'; e q ) (106) 
xx 
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and treat the <5-function term < x y\VV\4> > separately. Thus fP<p will be decomposed into W(j> + rV(f>. One obtains 

<xy\VV\cj>> = yd' M ^^) ^ 6 ^(-Y°V(x) 

ZA LM 

]T i Xo y/2\ + lC(l \ L;m Q ,0) 



J2V2kTT(-) k A^ loXo (x,y) [l + (- 



\A +A+fci 

v -r J-v, ; ^ZA.ioAoV^: L- 1 - ^ ^ J J 



/: 



1 3 1 

dt P k (t) Uno (\-x + y\)j Xo (q \-x- ~y\), (107) 



where A is given by (see Appendix B) 



dx dy y; >x L ' M (x, y)y k ° k (^ $)3&° (~* + £ |* - 



$L,L>SM,m A U L l \ (x,y)- (108) 



Now the projection onto the orthogonal states 3^ >' M (i, 2/) defines T^ 1 LM (x,y) 

^O' 1 »^ — / A~ J„~, -w* L > M (~ „-,\tO,1/ 



T^ LM (x, y)^Jdx dy y^ M (x, y)T^(x, y) (109) 
in terms of which we obtain the two pieces of the driving term T 0,1 

TiAm^v) = ^(h 3 ((-) lo + (-) l )SAL mo (110) 



7T y 3' 

dqq 2 Jx(qy) / dy'y' 2 / dy"y" 2 j x (qy') 



k A 



h aj \ (v'> y")C(l X L, to , 0) 

1 d*P fe (t)n(x, \±f + ly>\ ; e q )u no (\^f + 

V2S ALmo (-) lo V(x)J2^o 



k 

/ 1 i ^W^o(l^ + yl)jA (gol^-iy|)(i + (-) Ao+A+fe ) 

Here and in the following we use the notation s = 2s + 1. This symbol s should of course not be confused with the 
unit vector. 

n Before adressing the second part of the driving term in Ea H41|) or i|76[l connected to the 3- atom force we regard 
the part of the kernel without 3-atom force, < xy\tGoVT > . Using similar steps and projecting onto states of total 
angular momentum L leads to 

Kfcl (xy) = J dxd V yt x M *{xy) <xy\tG VT> (111) 

pOQ pOO pOQ 

= 4 / dqq 2 j x (qy) / dx" x" 2 / dx'x' 2 Ti(x, x' ; e q )gi(x' , x" ; e q ) 
Jo Jo Jo 

dy"y" JX (qy") £ £ Vk(-) 1 '+ k A^, x ,(x" y") 

dtp k {t)T vx , L <\\ar + f\, \- A x" - + (-f +x ) 
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where we have defined the quantities XJal in terms of the expansion 

T(ccy)= J2 y[{ M (xy)TiXL(xy) (112) 



IXLM 

of T(xy) given in Eas H41Jl and (|77(l . We see that the unknown amplitudes Ti\L{xy) occur under the integral with 
shifted arguments. 

We are left with two more terms in Eas. if4*T|) . an d ( T77|l . the second part of the driving term and another kernel 

piece, both containing the part v} 1 ' of the 3-atom force. Let us regard the simplest term < xy\V^\(j)i > explicitely 
We expand the plane wave and introduce states of total angular momentum, with the result 

<^yW|^> = ^0^Y.' X °&xM^)nn{x)V^{xyx-v) (113) 

Xo 

^C(l \ L,m ,0)y^(,xy) 

L 

Next we expand the 3-body force into Legendre polynomials 

V<U(x, y,x-y)=J2 ^K-) k ytl(xy)v k (114) 

k 

with 

v k (x,y)= J ^dtP k (t)V {1 \x,y,x-y) (115) 

and combine the angular dependent terms as 

1 



= ^kl \ (-) Xo+L+k (116) 

where the terms in curly bracket are 6-j symbols;^. Finally projecting E dll3() onto yl\ M (xy) one obtains 
J dx Jdy y l \ M *(m<^V (1) \<t>>=^^^i Xo Mjx (q y)u no (x) 



(117) 



Ao 

C(kl l, 00)C(fcA o A, 00)C(Z o A o L, m , 0) 



k 

/ dtPk(t)V l >(x,y,t) 

The other pieces of that second part of the driving term can be worked out similarily and we obtain altogether 

Tew = r *v r nisi 

jo jo xx yy 

pOO pOO 

+ ~ dqq 2 j x (qy)jx(qy') / dx" x" 2 Ti(xx" ; e q )gi(x"x'; e q )} 
n Jo Jo 

-=J2 lX °\ X oJ2" kC( - loX ° L > m °V / dtP k (t)V^(x'y't) 
[^jx (qoy>n (x , )\flX(-) Xo + L+l ^ ^ * \ |c(fcZ o /,00)C(fcAoA,00) 

(-)" 



^i_^C(fcfc' M ,00) 2 A^, io A (^y)(l + (-) An+A ' +A1+fc ')] 
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The remaining parts of the kernel in Ea H41|) or (|77|l can be worked out similarily and we split it into two pieces: 

K$£\xy) = JdxJ dyyt x M *(xy) < xy\V^G Q (l + V)T > (119) 

= -J2H-) k dtV^(xy,t)J2 dqq 2 j^(qy) 
n k J- 1 Jo 



in 

oo 



and 



dx'x' 2 / dy'y' j^iqy^g^ixx'ie^lT^^Lix'y') 
Jo 

+ J2(2n)Vk'(-fJ2 f dtP kl {t)T llVL {\\x' + $\^x'-\$\) 

((-)'' + {-) l ' +x ' + n^,vx>{x'y')\ 

K^ixy) = fdx f dyyff*(xy) < xp\hG V^G (l + V)T > (120) 



dqq 2 j\(qy) / dy'y' 2 j x (qy') / dx'x' 2 dx 



T l (xx";e q )g l (x"x';e g )K^\x'y') 

Then we end up with the one integral equation for Ti\i,(xy) in the form 

T lXL (xy) = T^ L (xy) + Tf& (xy) + (xy) + (xy) + K% 2) (xy) (121) 

where the T 0,1 are the driving terms and the if's incorporate the T-amplitudes. 

These are coupled equations among the T's, since Eqs ()lllll . ()119|l and (|120|l contain sums over several T's. Once 
the T-amplitudes in case of distinguishable particles are obtained we can use the Eqs. I|26|) . 127|l . (|28|) to determine 
the elastic and arrangement amplitudes as well as the complete break-up amplitude 120(1 . In the case of identical 
atoms one can not distinguish, of course, between elastic and arrangement amplitudes. Let us first regard the elastic 
amplitude given in Eg 14311 . 

For < cj)'\VT\(f) > one expands the plane wave, introduces the expression (|112f> for T, total angular momentum 
states and finally expands u n /(||y' -I- §j7|)T)AL(|§ y 1 + %y\>y') m t° Legendre polynomials. Then the same expression 
(|105|l for the remaining angular integrations appears and we obtain 

<0'|7>T|0> = 2V^) 3 E(-0 V *v mo -m'(<z') (122) 

6 A' 
poo poo 

/ dyy 2 d y 'y' 2 J2C(l'\'L,m'm -m') VVif-f 
Jo Jo L k 

y E h t$,ixM)((-) 1 + (-)'') 

;a 

^ dtP k (t)u n ,(\jy' + \y\)T lXL {\^ + ^y\,y') 

Here m' is the magnetic quantum number of the final 2-atom bound state with orbital angular momentum I' and h"'. 
has been defined previously in (|105|l . 
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Very similar steps applied to < (p'\VG 1 \4> > lead to 

<0'|PG o - 1 |0> = A(!)3£ H) A'y A , TO0 _ ro ,(<f) (123) 
>/5F 3 ^ 



rfyy 2 / dy' y ' 2 Y, iXo h-{<i'y)hM)\lH{-) 10 + (-)'') 

J] C(1'A'L, m'm - m')C(Z A £, mo, 0) ^ V^(-)*^, IA (tf, 2/') 

L k 

dtP k {t)u n ,(\^y> + + ^IKodfif + ^) 

There are two more terms in (|43|) which are worked out by similar steps as for the kernel parts. The resulting 
expressions are 

f \l poo roo 

<4/\V^(l + V)\<f>> = L>J2(~ifY x , m ^ m ,(q') / dxx 2 / dyy 2 u n ,(x)j x >(q'y) (124) 

V ^ y JO JO 

i^'v^oYl C ( l o x oL, m 0)C(l'X'L, m'm Q - m!) 

\ L 

J2'k f dtP k (t)V (1 \xyt) 
k 

[^\fi{-f V^onMoy)u no {x)(-) L { l Q x h \ }c(feZoZ',00)C(feAoA',00) 



J2k'(-f f dtP kf (t)u no (\^x + y\)j Xo ( qo (\lx- \y\)) 
k> 1 

J2 -^C(kk^,00) 2 Af y>loXo (xy)(l + (_)Ao+fe+V )] 



< (f>'\ (l+V)G T\(j)>= I dxx 2 I dyy 2 (125) 



W 3 / 2 Jo Jo 
u nl {x)J2{~i) x 'jx l (q'y)Y x , m( ^ ml (q , )J2C(l , X'L,m , m ) - m!) 

A' L 

/l />oo 
dtP k (t)V^(xy,t)Y] / dqq 2 j^(qy) 

E \/A^(-) M2+L+ " +i { Z 2 £ t } C(fc/x 1 Z,00)C(fc M2 A,00) 

/At ** ' 

poo 

dx'x' 2 / dy'y' 2 j^{<iy')g»A xx ') e q)\TiJ.i^L{x'y') 
o Jo 



£(27r)\/^R fe 'E /' dtP k ,{t)T vx , L {\\x> 

k' l'\> 



2/ I) \-^ x ~ -jV 



((-)'' + R Z ' +A ' +M2 )A^ a , z ,vW)] 
Finally the complete break-up amplitude < 0o|(l + "P)T > according to Eq. I|42l) can be based on its first term 

T(pq) = <fa\T> = <pq\T> (126) 
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Expanding the plane waves one readily obtains 



T(pq) = lj2ya m °m)(-i) l+X (127) 

IXL 

/>0O 

dxx 2 I dyy 2 ji(px)jx(qy)Ti XL (xy) 
o Jo 

The remaining two terms in (|42|) simply require to change the values of p and q as given in one example inf I85|l . 

This concludes the display of the formalism. The cases for two identical atoms out of the three or three distin- 
guishable atoms are somewhat more tedious to be worked out, but straightforward and are left to the interested 
practioner. 



V. THE S-WAVE EXPRESSIONS 

For very low energy of the incoming atom and the molecule in its ground state one can restrict the treatment to 
setting all orbital angular momenta to zero. Then the final expressions simplify and allow an overview of the structure 
of the integral equation to be solved and the quadrature expression for elastic scattering. Again we display only the 
case of identical atoms. In this case the quantities A and h given in the Appendices both simplify to and the 
integral equation l|121|) composed of the parts (jllUfl . (|118|) . (|lll|) . l|119|) and (|12U[) now reads 

T 00 o(x,y) = ^(-) 3 / dqq 2 j0 (qy) \ dy'y' 2 j (qy') \ dy"y" 2 j0 (q y") (128) 
l 



o Jo Jo 

4 2 2 4 



^ dt r (x, \^y" + -y<\;e q )u no (\-y" + ^\) 

If 1 ! 3 1 

+ -^~ V ( X ) dtu no (\-x + y\)j (qo\-:X - -y\) 



V2ir 



dx'x 12 1 w ^ -*')%-</) 



o Jo ' xx' yy' 



dqq jo(qy)jo(qy') / dx"x" t (xx"; e q )g a (x"x'; e g )] 
o Jo 



-L £ dtVM(x'y>t)[^=j ( qo y>)u no (x') + J* dtu no {\\x' + ^|)j ( 9o ||^ - \#\)^\ 

dqq 2 j (qy) / dx" x" 2 I dx'x' 2 T (x, x'; e q )g (x' , x"; e q ) 
it Jo Jo 

dy"y" 2 jo{qy") f dtT mQ {\\x" + y"\,\-/' \f\) 

^ rl poo noo noo 

+ - / dtV^(xy,t) / dqq 2 j {qy) / dy'y' 2 j (qy') / dx' x' 2 g {xx' ; e q ) 



" J-l JO Jo 

/! i 3 1 

^dtToooQ-x 1 +y'\,\-x' - -y 

2 

7T 



dqq 2 J (qy) / dy'y' 2 j (qy') / dx'x' 2 _O ac dx"x" 2 T (xx"; e q )g (x"x'; e q ) 
o Jo Jo 



i pi poo poo poo 

- dtVW(x'y',t) dq'q' 2 Jo (q'y') dx'" x'" 2 dy">y"' 2 j (qy"')g Q (x'x">; e q .) 
7f ./ — i ,/n ./n ,/n 



-1 JO JO Jo 

\rr { III lll\ i / lirri (\ ^ -'III . -Ill I I ^ ->lll ^ 

[Tooo(x y )+ dtT oo{\-x +y \,\-x - -y 
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The elastic amplitude given by the expressions (|119fl , ( 1120(1 , i|121|) , i|122|) also simplifies and one obtains 

i />OG />OC 

< <t>'W\<t> > = TTT / d VV 2 / ^y 2 iA'(9'2/bo(9oy') (129) 



2^ 2 7o 

l 



o 



^ dtu« o (|-0' + l^Wd^y 1 + ^KoG^ + -y\) 

fO a r°° r°° /*! AO OA 

^) 3 j o dyy 2 j dy'y' 2 J dtu no (\-y> + -y\)T 000 (\-tf + -y\,y') 



1 1 



dxx 2 / dyy 2 u n >(x)j (q'y) / dtV^^xyt) 



l 



l 



2 (2tt)3/2 7 7 

/ 1 1 3 1 

dtu no (\-x + y\)j (q (\-S - -y\)] 

-7=^2 dxx2 / dyy 2 u n >{x)jo(q'y)—= / (ftF (1) (:r2/,t) / dqq 2 j (qy) 

dxx 2 I dy'y' 2 j (qy')g (xx';e q ) 
o Jo 

Pooo(*V) + I' *Tooo(|^ + y'l 1^' - ^1)] 
At this low energy complete break up will not be possible. 

VI. COMPUTATIONAL CONSIDERATIONS 

The numerical solution of the integral equations for the T-amplitudes requires a large amount of computer time, 
depending on the number of angular momenta that enter in the partial wave expansion, and also on the number 
of mesh points that are needed for the discretization of the integration kernels, or equivalently, on the number of 
mesh points needed for the discretization of the coordinates x and y. We expect to minimize the number of required 
mesh points by using a recently developed spectral integral equation method (S-IEM) Q for solving the two-body 
Lippmann-Schwinger integral equation in configuration space that is very economical in the number of mesh points 
required for a given accuracy. This feature has now been demonstrated in several applications that are summarized in 
. For example , for the case of the binding energy of the He-He dimer which is very small, the calculation of the 
bound state wave function has to be carried out to large distances. For a distance of 3000 a.u. the S-IEM required 
only 200 mesh points to obtain an accuracy of three significant figures for the binding energy, and with 320 mesh 
points the accuracy increased to six significant figures Q. The method consists in dividing the radial interval into 
partitions, each partition receiving a fixed number of mesh points, expanding the unknown wave function into a series 
of Chebyshev polynomials in each partition, and then solving for the coefficients of the expansion. Such expansions 
are by themselves very efficient, and their accuracy properties are known. In addition, the size of the partitions is 
made automatically small in the region where the wave function changes rapidly, and large where it changes slowly, a 
feature that further contributes to the economy. An additional advantage of the spectral method is that conventional 
interpolation methods can be avoided. In the Faddeev scheme the Jacobi coordinates in one arrangement have to 
be translated into the corresponding coordinates in another arrangement. These translated points do not fall onto 
the predetermined mesh of points in another arrangement, hence interpolations are usually required with methods 
using fixed mesh points. The spectral method avoids this problem, since the Chebyshev polynomials, being analytic 
functions, can be evaluated at any prescribed positions, and hence provide the necessary translations. 

The two-body r - matrices play a large role in the calculation of the driving terms (II 10B ( I118|) and the integration 
kernels <|lllO ( I119|l f ll20|) . Investigations in progress are showing that the S-IEM can also be applied to the evaluation 
of the two- variable r- matrices, and the partition structure of the radial intervals of the latter will then subsequently 
determine the partition structure of the integration kernels. The number of mesh points for each of the two variables 
in the r -matrix is expected to be of the order of 50. The size of the matrices that represent the integration kernels 
will have to be estimated, and depending on the outcome, iterative procedures may be required. 
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VII. SUMMARY AND CONCLUSIONS 



The purpose of this paper is to lay out the theoretical formulation of the solution of the three-body Faddeev integral 
equations in configuration space, for the purpose of applying them to atomic physics situations. It is assumed that 
the Born-Oppenheimer approximation is valid and the three atoms are moving on one potential surface. The case 
of a conical intersection is not treated but it appears conceivable that also such a case could be handled by using 
coupled Faddeev equations, similar to the treatment of A - excitations in a three-nucleon problem. We have chosen 
a formulation of the Faddeev equations which turned out to be extremely useful in nuclear physics 0. Instead of 
working with wave function components we introduce break-up amplitudes T,;, which vanish at large distances and 
fulfill a set of three coupled equations. In case of identical atoms only one equation is needed. That set of equations 
is then displayed in a configuration space vector representation and furtheron decomposed into partial waves. For 
the sake of clarity this set of coupled equations is also presented assuming s- waves only. Once the T-amplitudes 
are determined all physical matrix elements for elastic and inelastic atom -diatom scattering, for the arrangement 
processes and the complete 3-atom break- up process are obtained by simple quadrature. The relevant expressions 
are given in detailed form. 

One reason for working in configuration space is that the atom-atom two- and three-body potential surfaces are 
given in that space. One objection commonly raised against performing three - body calculations in configuration 
space is that one needs many mesh points for the numerical implementation, because the wave functions have to be 
calculated out to large distances in order to impose the asymptotic boundary conditions. We overcome this objection 
for several reasons. One is, as already said, that we do not calculate wave functions, but rather T-amplitudes. 
These functions are essentially the product of wave functions times potentials, and hence decay fast with distance in 
contrast to wave functions, which oscillate at infinity. Also, the asymptotic boundary conditions of the underlying 
wave functions are automatically included, because the Faddeev integral equations contain the Greens functions which 
lead to the appropriate asymptotic behavior. Another objection against working in configuration space is that wave 
functions, and so also T-amplitudes, can have a strongly oscillatory behavior, and hence many mesh points may be 
required. We expect to overcome this objection by using a recently developed spectral method (S-IEM) for solving 
the two-body Lippmann-Schwinger integral equation that is very economical in the number of mesh points required 
for a given accuracy The economy in mesh points is crucial when solving the Faddeev integral or differential 
equations, because the numerical complexity increases like the cube of the dimension of the final matrix, and the size 
of the dimension is proportional to the number of mesh points. 

In summary in this paper we lay out in detail the partial wave expanded Faddeev equations for the T-amplitudes 
in configuration space and in integral form. These equations, although complicated, are not much more complicated 
than the corresponding equations in momentum space, which have been successfully solved in the realm of nuclear 
physics. In a future study we intend to solve the equations for a simple test case in order to study the numerical 
feasibilty of the method. 
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APPENDIX A: THE EXPRESSION^, v (y, y') FROM EQ (tl05l) 

The evaluation of h can be done using standard angular momentum algebra [{| So we do not provide the 

detailed steps but mention only some useful formula, which are needed and possibly not so commonly in use. 

The dependence on the two vector directions in the shifted arguments for the spherical harmonics can be seperated 
using 



Further in the course of recoupling the following formula is needed 
The resulting expression for h is 



/ (2/ 1 + l)(2/ 2 + l) 
4tt(2Z + 1) 



C(hl 2 l,00)Y lm (a) 



ht&vM) = jdyj dy'y;x +^my^v')yiS:^ +^v') 
^ E E v u+ y l2+l '^r + H^ 



h+l 2 =l l[+l' 2 =l' 

1 (2Z + 1)! / 

(2Ii + l)!(2I 2 + l)!V (2Zi + 1)1(2/^ + 1)! 



(21' + 1)1 



E 



IhX^J l%\> \J kl[ 

h h I 
A L f 



C^A/, 00)^| Jj \ *' \C(l' 2 ,X'f, 00) 



E(-) 7l+fc {/' I' ~ x }c(fc/ l 7i,00) 2 C( 7l Zi/,00)C(7iA,00) 

71 1 



APPENDIX B: THE EXPRESSION A^ lly (xy) FROM EQ ltT08l) 



A a L i>\> i x v) = J dxdyy* x L ' M (x, y^fix, y)y^(^S + y, - ^y) 

= ^ E E (~)^(f)^-) A ^+*v^ 



(2Z' + 1)! / (2A' + 1)! 



V (2iJ)!(2£! V (2A' 1 )!(2A^! 

A; \\ A' lc(ZiA / 1 //i,00)C(4A^2,00) 

^ V7I(-) /l yMcQehm, ooo) VmM-r +i+71 (-)' 

h 71 

C( 7 iMiZ,00)(7(7iM2A,00) 
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